Pion correlations in Nuclear Matter 



P.K. Panda,/ S. Sarangi,,^ and J. da Providencia^ 

^Centra de Fisica Computacional, Departamento de Fisica, 
Universidade de Coimbra, P-3004 - 516 Coimhra, Portugal 
^ICFAI Institute of Science & Technology, Bhubaneswar-751 010, India 

Abstract 

The saturation properties of the nuclear matter taking pion correlations into account is studied. 
We construct a Bogoliubov transformations for the pion pair operators and calculate the energy 
associated with the pion pairs. The pion dispersion relation is investigated. We next study the 
correlation energy due to one pion exchange in nuclear matter and neutron matter at random 
phase approximation using the generator coordinate method. The techniques of the charged pion 
correlations are discussed in the neutron matter calculations. We observe that there is no sign of 
the pion condensation in this model. 
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I. INTRODUCTION 



The understanding of the nuclear force at a microscopic level is an important problem 
since it shall be the basis for the nuclear matter and finite nuclei calculations. The inter- 
actions of nucleons which may arise as a residual interaction due to their substructure of 
quarks and gluons is technically not solvable. The alternative approach is to tackle the 
problem through meson interactions. 

In recent years relativistic mean field theory (RMF) has been quite successful in 
describing nuclear matter and finite nuclei properties. The NN dynamics arises in this 
model from the exchange of a Lorentz scalar isoscalar meson, cr, which provides the mid 
range attraction and an isoscalar vector meson, cj, which provides the repulsion. Also, in 
this model the inclusion of the p- meson takes care of the neutron-proton asymmetry. With 
a small number of parameters, the RMF model reproduces the nuclear matter saturation 
and describes the bulk and the single particle properties for the finite nuclei reasonably well. 

Despite the success of the RMF model, several open questions still remain unanswered. 
Firstly, the meson fields are classical and secondly the attractive part of the nuclear force is 
mediated through the hypothetical a meson which could be an effect of multi-pion exchanges 
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3|, U, 15|, la, LZl]. The cr can not be interpreted as the representation of a physical particle, 
since such a particle or resonant state remains still to be confirmed in the experiment. Thus 
on aesthetic as well as phenomenological ground, alternative approaches will add to our 
understandings. 

The original Walecka model at the Hartree approximation does not contain a dynamical 
description of the pion fields. However the importance of the pions in NN dynamics can 
not be ignored. Realizing the essential role of pion in the description of the nuclear medium 
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an alternative approach for the nuclear matter deuteron [5| and for ^He p] has been 
developed. In this method, it has been studied the description of nuclear matter using pion 
pairs through a squeezed coherent state type of construction 4,[5, 6, 7|. This simulates the 
effects of the a-meson and is a very natural quantum mechanical formalism for the classical 
fields. 

The generator coordinate method (GCM) is a technique of great physical appeal which has 
been developed jsl to describe collective oscillations in nuclei. Besides being extensively used 
in nuclear structure physics, it often finds application in various other branches of physics 



9j . In this work we investigate the pion condensation problem in nuclear matter and neutron 
matter using the generator coordinate method. In a simplified model, the problem had been 



studied in Ref. [ly] where a coherent state description for the pions was used. Using the 



gaussian overlap and harmonic approximation, the Hamiltonian may be diagonalized by a 



randam phase approximation (RPA) like canonical transformation ll|] . In this methodology, 
one can go beyond the coherent state description. Again the calculation can be carried out 
exactly, without further approximation. The present analysis is an extension of the mean 
field approach of Walecka where classical fields are replaced by quantum coherent states for 
the pion pairs. This report has also an advantage that the one pion exchange correlation 
contributions are considered at RPA level using similar Bogoliubov transformations. In the 
present model, we have observed no sign of the pion condensations in the RPA modes. 

The outline of the paper follows. In section 2, we derive a pion nucleon Hamiltonian in 
a non-relativistic limit. We then construct a Bogoliubov transformations for the pion pair 
operators and calculate the energy associated with the pion pairs. In section 3, we calculate 
the correlation energies due to one pion exchange in nuclear matter and neutron matter at 
random phase approximation (RPA). Section 4 consists of the discussions of the saturation 
properties of nuclear matter, the pion dispersion relation in the medium and a concluding 
remarks. 

II. FORMALISM 

A. Non-relativistic Hamiltonian 

The Lagrangian for the pion nucleon system is taken as 

£ = ^ (^^^5^ -M + G^^ct)) (d^^.d^^, - m\iip,) , (1) 
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is the doublet of the nucleon field with mass M, yjj's are pion fields 

and = Ti(fi represents the off-mass shell isospin triplet pion field with mass m. G is the 
pion-nucleon coupling constant. Repeated indices indicate summation. 
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Prom the above Lagrangian, the equation of motions are 



{E - M)'^i -{a-p + iG(f))7pn = , 



(2) 



{E + M)^H - (a ■ p - iG(t))^i = , 



(3) 



where, E — i{d/dt) and p — i{d/dx). Ehminating the small component ■0// from equation 
(2) and (3) we have 



[E^ -M^) - {E + M){a-p + iG4>){E + M)-\a-p-iG(j))] ipi ^ . 



(4) 



Equation (4) can be rewritten as 



E'-M' -p' + iG[{a ■p),(t)]-G'(t)-(t) ^i = Q 



(5) 



From equation (5), we can immediately identify the effective Hamiltonian for the nucleons 



as 



(6) 



where the single particle nucleon energy operator e^, is given by = (M^ — V^)^/^. In the 
non-relativistic assumption, we have to replace by M, when in a denominator and by 
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M + when not in a denominator. Now the effective Hamiltonian becomes 



(7) 



where the free nucleon part H^(x) is given by 



2M 



(8) 



the ttN interaction Hamiltonian is provided by 



iG 
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(9) 



and the free meson part 7ij\/(x) is defined as 



(10) 



We expand the pion field operator ipi{x.) in terms of the creation and annihilation operators 
of off-mass shell pions satisfying equal time algebra as 
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[ai(x)"^ + ai(x)), 



V9»(x) = i^l y (a,(x)T - a,(x)) , (11) 



with energy = {rn? — V^)^/^. 



B. Correlation energy associated with two pions and Bogoliubov Transformation 

The quadratic terms in the pion field in eq. ([9]) provide a isoscalar scalar interaction of 
nucleons and thus would simulate the effects of cr-mesons of the Walecka model. 
A pion-pair creation operator given as 

i^^ = ^E/k4.«lk., (12) 

is then constructed with the creation and annihilation operators in momentum space and 
the ansatz function /(k). We then define the unitary transformation U as 

U = e(^'-^) (13) 



and note that t/, operating on vacuum, creates an arbitrarily large number of scalar isospin 



singlet pairs of pions corresponding to squeezed coherent states. We will show that this is 
the appropriate transformation to diagonalize the pion part of the Hamiltonian. The "pion 
dressing" of nuclear matter is then introduced through the state 



1^) = t/|0) =e(^'-^)|0). (14) 

We obtain 

Ski = a^iU = (cosh/k) aki + (sinh/k) al^.^, (15) 

which is a Bogoliubov transformation. Here f/ is a unitary and hermitian operator. The 
psedo-pions ciki are the results of the unitary transformation. It can also be easily checked 
that the operator a^i satisfies the standard bosonic commutation relations: 

[Ski, al'j] = ^jj4,k', [aii, al,j] = [a^i, Ok'j] = . (16) 

and also 

aki|^) = (17) 

The reverse transformation: 

Oki = (cosh /k) ttki - (sinh /k) alki = a^i - l/k a^u ■ (18) 
In momentum space the effective Hamiltonian may be re-written as 

pqjaa'rjrj' ^JVl^UJqV 
pqjQj; q 

Here p, a and ?7 are respectively, the momentum, spin and iso-spin quantum numbers of the 
nucleon and q, j are the momentum and isospin labels of the pion. p = |p|, q = |q|. cj,^^ is 
the creation operator for nucleon with momentum p, spin a and isospin rj. The contribution 



of the quadratic term in the pion field coming from the above Hamiltonian as 



+ E 



^q«qj"qj + (,"q,i"-qj + «qJ«-qJ + ^Sj^q.J j 
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where = (wq + ^ j = c^q + (?' with g - 



Gin] - ^ n' ..rHh n'- G1p_ 

Now the equation of motion for the pions becomes 



This gives 



q '^qj ~l~ 9 -^q '^qj ^qUq ^qj 9 Vq ^qj ^q{^q(^qj ~l" Vq^—qj)- 



The characteristic equation is 



~2 / ^ I /2 



which gives 



'^q = - g''^ with 



2Un 



Now 



i^2^ = E'^q«qj«qj + o E('^q " ^q) • 



^q"'qj"'q,> ■ 9 
qj ^ q 



We now have to include a term which corresponds to a phenomenological repulsion 



between the pions of a "pair" in the above Hamihonian if27r and is given by 



H^ = AY: e^^'^'a^aq, (26) 

where the two parameters A and i?^ correspond to the strength and length scale, respectively, 
of the repulsion and will be determined phenomenologically. This term amounts to imposing 
a cut off on the momentum q which accounts to the fact that momenta larger than kf are 
not dynamically meaningful. With this repulsion term, we now have 



< = ( ^q + ^ e^^'^^ + ^ 1 = + ^ ^""^'^ + 3' (27) 



with 



After transformation the Hamiltonian in equation f[T9l) becomes 
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5^ / '^p+qjOryCp^a'^' (zcr • q)aa'Tj(aqj + Ct-^qj) (2^) 



pqjaa'riri' 



where 



_ g (Xq - 1/q) 
^?q- 2M • ^^"^ 



In the next section, we will consider RPA and calculate the correlation energy associated 
with one pion exchange for nuclear matter and neutron matter. 

III. CORRELATION ENERGY ASSOCIATED WITH ONE PION EXCHANGE 
A. nuclear matter 

We consider a Slater determinant of plane waves 

1"^)= n (31) 

aV,\p\<PF 



with |0) is the absolute vacuum, is the Fermi momentum and 



Cp,a,.|0) = . (32) 

Excitations with momentum transfer q are coupled to excitations with momentum transfer 
— q. Thus, the wave function |^') which describes such excitations should read 

|*)=exp5|$), (33) 



where 



a,a' ,ri,Ti'p € f2q 
a,a' ,Ti,ri'p G f2-q 

= U^K. + ^-.^-q. • (34) 
In the above, J\fq is the normalization factor insuring 

Yl \{a,v\a-cirj\a',ri')\'^4M, E = 1 (35) 

aa'rjrj' ,p G f2q p G fiq 

and the domain is defined by |p + q| > Pfi |p| < \pf- Only positive energy states are 
occupied. Here, ja,?]) denotes the spin iso-spin eigenstate, a^\a,vi) = a\a, rj) and Ts\a,ri) = 
7]\a,r]). The determination of the A/"g is now very simple. All we need is the volume of the 
intersection of 2 spheres of radius pp, theirs centers being a distance q apart. With this 
assumptions, we have 



- (2^'"^ - 12 
The transformed unperturbed Hamiltonian becomes 



K' = 7^.'^-l'{A--] ■ (36) 



Ho=J2^P 4>,ar,Cp,ar, + E ^Aj<i + ^ ' ^ = o E('^9 " ^q) " (37) 

P,otV qj 9 



The pion nucleon coupling reads, 



Hint = - ^ {ar]\i{a ■ q)Tj\a'r]')^^ cJ,+q,a'?Cp,aV «j + «-qj) ■ (38) 



In order to proceed, it is convenient to bosonize the Hamihonian H, restricted to the sub- 
space 5". This is done by the replacement Sgj Bg j satisfy boson commutation relations. 
The bosonized Hamiltonian reads as 

qj 

The parameters of Hb, namely Eg and Qg are fixed by expectation values such as 
{^\B^HBIj\^) = Ehf + Eg, {^a^HB]^\^) = Qg, 

and 

m~a^B_^H\^) = Qg. (40) 
The above Hamiltonian is diagonalized by a Bogoliubov transformation of the type 

= + + yt^B_^ + yt^~a.g,, n = 1, 2 (41) 

which leads to excitation energies and the correlation energy. 
In the above 

Eg = AMgq^ J2 TT^iiP^^ -P)- 



1 4gp|, 



2M3(p|n-gVl2) ' 

and 

The eigen frequencies are 



\2UJgV 



(42) 



^gg . (43) 



^^t^ = + (^l ± /(£^ - (^If + l&Sgu'^g . (44) 



The correlation energy becomes 



B. neutron matter 

The correlated Fermion wave function may be written, r = — 1 for neutron and r = 1 for 
proton 

|*) = exp5|$), n 4a,-i|0), (46) 

a,|p|<PF 

where the correlation operator reads 

S = U^Bl, + U_^Bl^^o + V^Bl^ + , (47) 

with the quasi boson operators 

K,0 ^ \fK E (">-l|(c^-q)^o|a',-l) 4+q_«^_lCp,a',-l (48) 
a,a',p € Qq 

^ E («>l|('^-q)'^+l«'>-l) 4+q,MCp,a',-l , (49) 

a,a',|p| < Pf 

where Tq = T3 and t+ = (ri + ir2)/2. Here, A/'g and J\f'g are normalization factors insuring 

Af, E |(a,-l|a-qr,|a',-l)|2 = 2A4 E = 1 ' (50) 
AA', E |(a,l|a-qr,|a',-l)r = 2Ar', E = 1 • (51) 

a,aMp| <PF P <PF 

As earlier the domain fig is defined by |p + q| > Pf, |p| < \pf- The determination of the 
normalization, J\fg and jV^, is now very simple. To compute Afg, all we need is the volume 
of the intersection of two spheres of radius pp, theirs centers being a distance q apart. We 
found 



The kinetic energy for the particle-hole pairs with momentum q become 
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4 = 2A/;g^E^((p+^)'-p') 
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P<Pf 

The pion nucleon interaction becomes 

Hint = - X! {<^V\K^ ■ q)'ro|a'^')-^==cJ,+q,a,?Cp,a',7'(«q,0 + olq.o) 
pq,aa'r?T)' y ''^ 

- Yl {C^VW ■ q)^+l«V)^==4+q,a^Wr?'«+ + «-q,-) 
pq,aa'r;?j' y ''^ 

- {(^nW ■ q)^- l"V) P—A +ci,ar,Cp,a'v' «+ + «-q -) 



The effective bosonized Hamiltonian containing two pion exchange becomes 



Hb = Eq + ^1 £g5j o5q,0 + ^q^i+^q,+ + '^q X! "qj^qJ 

- E<5g(<0 + 5q,o)«0 + alq,o) 

q 

- 2Z + «-q -) + ^q,+ («-q - + 4,+)) ' 



where 



Q'a 



9q 



The RPA equations are easily obtained. For the modes with charge, 

[Hb, {XqBl+ + Cgai,+ + Vqa-q,-)] 

= Bl^+{Xqe'q + CQ'q - TjqQ'q + 4,+ (^q<3g + Cq^q) + ^-^A^qQ'q " Vq^q) 



The characteristic is a cubic and becomes 



g; (Dq - 



-(Dq — Qg 



= -Q3 + £^q2 ^ ~ ^ 2Q'^cI;q - s;^;^ = (59) 



/2~ /~2 



where the solutions are 



e 

3 + 



21/3 _ 3~2j 



1/3 



(60) 
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(1 + iV^) - 3ou^) 



3 X 22/3 (^-2e'^ - 54Q'^Co + ISe'cD^ + ^4 (-e'^ - 36^2)"'^ + (^-2e'^ - 54Q'^(Jj + ISe'Cj^f^ 



1/3 



6 X 2V3 



^(3) _ £: (l-zV3)(-£-^-3^2) 

1 ~ 3 / I 3 2\^/^ 

3 X 22/3 -2£'3 _ 54g'2cD + 18£'cu2 + J4 (-£'2 - 3c^2y + (^_2£'3 _ 54g/2^ + 18£'(I;2) 

+ ^ ^"'21/3 (1 + ^^^) (^-2^'^ - 54g'^c:' + ISe'cD^ + ^4 (-e'^ _ 3c^2y'^ + (^_2£/3 _ 54g/2^ + is^', 



Similarly 

[//b, (FgS_q,+ + Cqa-^,+ + ^y^Sq -)] = ^q{YqB-^,+ + C9alq,+ + TygCq _) . (63) 



This leads to a similar equation with replaced by —^q, so that the eigenfrequencies occur 
in pairs, ±Qq. The correlation energy for the charged modes becomes 

ELrr = I Ed^i'^ I + l^f I + l^f ^ I - < - 2^.) • (64) 

^ q 

The eigen frequencies of uncharged modes are 

= ^ + ul ± ulY + IQSqUqQl . (65) 

The correlation energy becomes 

Ecorr = \^{^^^^ +^[-^ -eq~Cjq) . (66) 

which are given in earlier section. 



IV. RESULTS AND DISCUSSION 



We first proceed to describe the binding energy for nuclear matter. We obtain the free 
nucleon kinetic energy density 

hj = {^\Tr[p^nd^)m = 1^ [M+-jlj . (67) 

In the above equation, spin degeneracy factor 7 = 4 (2) for nuclear matter (neutron matter) 
and , kf represents the Fermi momenta of the nucleons. The Fermi momenta kf and the 
nucleon densities are related by fc/ = (Qn'^p/'y)^ . It is well known that the short range 

TABLE I: Parameters of the model obtained self consistently at saturation density. 



a R-iY 
(MeV) (fm) (fm^) 
14.58 1.45 3.07 



interaction plays a crucial role in determining the saturation density which is mediated by 
the iso-scalar vector 00 mesons. Here we introduce the energy of repulsion by the simple 



form 



A.p^ (68) 



where the parameter is to be fixed using the saturation properties of nuclear matter as 
described in Ref. j?]. Thus we finally write down the binding energy per nucleon Eb of the 
symmetric nuclear matter (SNM): 



E 

Eb = — -M (69) 



where 



3 

Eo = hf + -J2^-i + hu.- (70) 
^ <? 

In the above, Eq is the energy density of nuclear matter or neutron matter without one pion 
correlation. The expression for Eq contains the three model parameters a, Rt^, and as 
introduced in the earlier section. These parameters are determined self-consistently through 
the saturation properties of nuclear matter at saturation density po = 0.15 fm^^ with and 
without the correlations. While pressure P vanishes at saturation density for symmetric 
nuclear matter, the values of binding energy per nucleon are chosen to be —16 MeV. In 
the numerical calculations, we have used the nucleon mass M = 940 MeV, the pion masses 
m = 140 MeV and the omega meson mass, = 783 MeV, and the ir — N coupling constant 
GV47r = 14.6. 




FIG. 1: Binding energy of symmetric nuclear matter (SNM) and pure neutron matter (PNM). The 
correlation is related to one-pion exchange 
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FIG. 2: The correlation energy from one-pion exchange in symmetric nuclear matter (SNM) and 
in pure neutron matter (PNM) 



We now discuss the results obtained in our calculation. We first construct a Bogoliubov 
transformation for the pion pairs operators and calculate the energy associated with it. 
We next calculate the correlation energy due to the one pion exchange in nuclear matter 
and neutron matter at RPA using the generator coordinate method. The binding energy 
per nucleon as a function of the density of the system is often refered as the nuclear 
equation of state (EOS). In figure 1, we present the EOS with and without correlation for 
the nuclear matter and for neutron matter. The correlation is related to one pion exchange. 
As expected, the binding energy for nuclear matter with and without correlation initially 
decreases with density and reaches a minimum at p/ pQ = l and then increases. 

In figure 2, we show the variation of the correlation energy, Ecorr as a function of density 
for symmentric nuclear matter (SNM) and for the pure neutron matter (PNM). The corre- 
lation energy initially decreases with density and then increses after the saturation density. 
The correlation energy for the neutron matter gives larger as compared to the nuclear matter 
at different densities. 

Dispersion relation of modes with the quantum numbers of the pions in nuclear medium 
is an interesting aspect. In figure 3, we plot the dispersion relations arisig from the two pion 
coherent states versus with the momentum at different densities. The increase of the pion 
dispersion relation for k/kf > 0.6 is probably an artifact of the repulsion term of equation 
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FIG. 3: The dispersion relation with the quantum number of the pions, a) for nuclear matter 
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FIG. 4: The dispersion relation of the RPA mode with the quantum number of the pions, r2_, at 
high densities (upper panel) and for low densities (lower panel) for nuclear matter. This corresponds 
to zero sound modes. 



(26). 

For nuclear matter, we observed two RPA modes with quantum number of the pions, 
fi-t. In figure 4, we have shown the dispersion relation for the RPA modes, versus 
momentum at different densities for nuclear matter. At p = po? ^- increases very slowly 
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FIG. 5: The dispersion relation of the RPA mode with the quantum number of the pions, r2+, at 
high densities (upper panel) and at low densities (liwer panel) for nuclear matter. It is showing an 
increase with density of the effective mass of the pions. 
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FIG. 6: The dispersion relation of the RPA mode with the quantum number of the pions, and 
I $73 1, for neutron matter 



with momentum. However at p = 4po, it increases fast. In the lower panel of the figure 4, 
we liave plotted VL_ for smaller densities. At small densities, f2_ increases monotonically 
with momentum k and corresponds to zero sound modes. 

In figure 5, we have shown the dispersion relation for the RPA modes, versus mo- 
mentum at different densities for nuclear matter. It is found the the magnitude of the f2+ is 



larger compared to In lower panel of the figure 5, we plotted for smaller densities, 
showing an increase with density of the effective mass of the pions. 

We next study the RPA modes for the neutron matter. There are three modes for the 
charge pions. In figure 6, we plot the RPA frequencies versus momentum k for neutron 
matter. The and \Q2\ are equal. In the upper pannel, we show versus momentum 
k. All the three RPA frequencies increase with density. In the lower panel, we plot IQaj 
versus momentum k which corresponds to zero sound modes. It is found that there is no 
sign of pion condensation in higher densities with the RPA modes. 

In conclusion, we have derived a pion nucleon Hamiltonian in a non-relativistic limit. 
We then have constructed a Bogoliubov transformations for the pion pair operators and 
calculate the energy associated with the pion pairs for the nuclear matter and neutron 
matter. This is an extension of the mean field approach of Walecka where the classical fields 
are replaced by the quantum coherent states for the pion pairs. We then calculated the 
correlation energies due to one pion exchange in nuclear matter and neutron matter at RPA 
using generator coordinate method. It is found that there is no sign of pion condensation 
in higher densities with the RPA modes. 
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